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Abstract 

The hypervelocity OB stars in the Milky Way Galaxy were ejected from the central regions some 10-100 million years ago. We argue 
that these stars, as well as many more abundant bound OB stars in the innermost few parsecs, were generated by the interactions of an 
AGN jet from the central black hole with a dense molecular cloud. Considerations of the associated energy and momentum injection 
have broader implications for the possible origin of the Fermi bubbles and for the enrichment of the intergalactic medium. 
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1. Introduction 

Hypervelocity stars (HVS) in our galaxy defy explanation. With 
velocities directed outward from the centre of the galaxy of be- 
tween 300 and 1000 km/s, these young OB stars cannot be ex- 
plained by ejection from binaries or via binary encounters with 
the central black hole. 

Most hypervelocity stars seem to be the relics of some- 
thing more exotic than binary ejections. The more massive the 
star, the higher the run-away fraction. Since a significant frac- 
tion are O stars of > 4OM , this means the event that gener- 
ated them was relatively recent. Binary star scattering by a cen- 
tral star cluster fails by at least two orders of magnitude to ac- 
count for their frequency dPerets & Subrfl2012l) . Binary ejection 
in supernovae seems unlikely given that their orbits are consis- 
tent with coming from the nucleus of the Milky Way Galaxy 
(MWG) and the travel times of hypervelocity B supergiants, 
around 3 - 4M , at 50-100 kpc from the Galactic centre (GC) 
and with velocities i 300km/s, are constrained to 60-200 Myr 
dBrown et al 1120121) . 

More exotic possibilities have been investigated. The Milky 
Way's central supermassive black hole might itself have a com- 
panion black hole, which could kick off the stars at high speed 
(lYu & Tremaind 12003b iBaumgardt etaD 120061: ISesana et all 
120081) . The main difficulties with this scenario are that the ve- 
locities could come out to be too high ( Ses ana et al ] l2008h . and 
that if dynamical friction with the stars in the dense Galactic 
centre is taken into account, a continuous supply of intermedi- 
ate mass black holes inspiralling into the GC must be invoked, 



because each one may o nly eject stars for a timescale of 10 6 yrs 
(IBaumgardt etalJl2006h . About 100 HVS (B stars of 3 - 1OM ) 
are generated in the halo over a typical propagation time of 100 
Myr, or at a mean ejection rate of 1/Myr ejected. We infer that 
the ejection most likely was spread over 10 6 to 10 8 yrs because 
of the spread in distances travelled and lifetime considerations. 

Recurrent AGN activity is an established piece of galaxy 
evolution, in contrast to recurrent accretion of intermediate mass 
black holes. Hence if there is no other convincing explanation, it 
is useful to estimate whether the mechanism of origin of hyper- 
velocity stars may be due to positive feedback on molecular gas 
induced by jet interactions generated by past explosions from 
our central supermassive black hole (SMBH). 

AGN jets have recently been shown via 3-D simulations to 
overp ressure clouds and in duce star formation in gas-rich central 
disks dGaibler et alJl201lh . Another recent study of AGN jets in- 
teracting with an cloudy interstellar medium demonstrates that 
10% or more of the jet energy effectively accelerates the gas 
clouds to escape vel ocity over a few tens of millions of years 
dWagneret al.N2012h . Neither of these studies includes cloud 
self-gravity, and one must await new studies before reaching any 
definitive conclusions, for example on cloud survival. However 
it is nevertheless useful to point out that there may be a local 
counterpart of jet-induced gas flows and induced star formation, 
for which the energetics and efficiencies seem to work out sur- 
prisingly well. 

Our own central black hole currently has a very low ac- 
cretion rate, but this may not have been the case in the past. 
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Indeed, jet activity has r ecently been claimed from radio data 
(lYusef-Zadeh et aLll2012l) . This leads the possibility that an eus 
(AGN)-like phenomenon may have repeatedly occurred in our 
own Galactic centre. Indeed, the recently discovered Fermi bub- 
bles in the central regions of the Galaxy may provide evidence 
for an outburst some 10 million years ago, provided reacceler- 
ation of energe t ic electrons occurred in the associated shock s 
dSu et all 120 lOt iGuo & Mathewsl 1201 It IZubovas et all 1201 ll) . 
Hypervelocity B stars suggest an event some 100 million years 
ago. 

Our starting point is that of momentum considerations. There 
is an intriguing coincidence between the momentum flow from 
the GC and that in hypervelocity stars that connects, we believe, 
to an outstanding problem in galaxy formation theory. The point 
here is that production of HVS is inefficient. Most of the stars 
formed by jet-induced overpressuring of GC gas clumps are in- 
evitably at relatively low velocity. These are tracers however of 
the same mechanism that generated the rarer HVS. Hence we 
first describe star formation in the inner few parsecs. 

2. Star formation in the Galactic centre 

The Milky Way Nucl ear Star Cluster has a radius of 5 pc and a 
mass of 2 - 3. 10 7 M o (lLaunhardt et al.l2002l) . It has several com- 
ponents, including the Sgr A* star cluster which contains ~ 200 
young stars within the central parsec. There is also disk struc- 
ture, with a main and an inclined disk at 0.8-12 arcsec (1 arc-sec 
= 0.04pc) of massive stars formed 6 Myr ago, as well as a more 
isotropic distribution of older (> 1 Gyr old) s tars containing 90 % 
of the stellar mass within the central parsec dPfuhl et alJIiOllh . 

The IMF of the WR, O, B stars in the central disks (in the 
inner most 4000 AU, 0.8-12 arc-sec) is top-heavy dBartko et al.l 
120101) . The WR and O stars formed coevally about 6 Myr ago. In 
addition there are late B (less massive) stars more isotropically 
distributed. Closer in there is an old star cluster with a standard 
IMF centred on Sgr A*. The B stars beyond 12 arc-sec also have 
a standard IMF. There could be a correction for less massive 
stars, possibly as large as 10, in the central arc-sec. If the central 
OB stars (within 1 arc sec) were also formed by the jet interac- 
tion 10 Myr ago, the HVS ejection efficiency is inferred to be 
of order 1-10%. 

Further out, the Arches cluster alone ha s a pre sent day kine- 
matic mass 1.5 x 1O 4 M dClarkson et all 1201 lb . It is 26 pc 
from the GC and only 2 Myrs old. The precursor molecular 
cloud mass of an Arches-like young massive star cluster is ob- 
served to be at a typic al scale of 3 pc and mass of 10 5 M o , see 
iLongmore et all (120121) . In our model, the jet-driven bubble ra- 
dius accelerates and compresses gas out to ~ 30 pc, where it can 
continue to induce star formation. 

At 10% star formation efficiency (SFE), we would expect 
1O 5 M in stars. In fact we see some 3.1O 4 M in low mass stars 
within 0. 1 pc by using the integr ated light at HST resolut ion of 
the central arc-second of the GC (lYusef-Zadeh et aU2012l) . 

2.1. A giant molecular cloud around the central SMBH 

At high redshift, the specific star formation rate is observed to be 
high. This is due either to an increase in the efficiency of star for- 
mation related perhaps to the increase in gas supply or to a new 
mode of star formation. The former interpretation, if coupled to 
the metallicity dependence of molecular gas formation, fails to 
simul taneously fit low and high redshift specific star f ormation 
rates dKhochfa r & Silk 201 h lKrumholz & Dekelll20T2l) . 



The occurrence of a new mode is plausibly related to the 
rapid rise in AGN activity, that parallels star formation rate his- 
tories, by z ~ 2. AGN interactions with gas-rich disks can pro- 
vide positive as well as negative feedback, with positive feed- 
back naturally augm enting the specific star formation rate, as in 
iGaibler et al.l (1201 ll) . The latter study shows how the pressure 
enhancements associated with a powerful jet can pressurise the 
gas-rich disk and thereby trigger and accelerate the star forma- 
tion rate. One consequence is ejection of gas clumps. We argue 
here that simple considerations of the momentum transfer and 
energetics support the case that stars forming in the ejected gas 
clumps can account for the hypervelocity stars, if the last such 
jet episode occurred some ten million years ago at the GC, and 
followed other similar episodes with an appropriate duty cycle 
determined by replenishment of the gas reservoir. 

Dense molecular gas is mapped at the GC. Some is infalling 
and some is undergoing tidal disruption, but a significant frac- 
tion end s up in a central circumnuclear gas disk at 1 .5-4 pc from 
the GC dLiu et alJl2012h . The mass of the inner giant molecular 
cloud (GMC) amounts to 1O 6 ~ 7 M . The residual gas mass de- 
pends on the star formation efficiency and on the assumed IMF. 
Gas replenishment will guarantee repeated episodes of AGN ac- 
tivity and star formation. If we assume that non-axisymmetric 
gravitational instabilities drive GMC formation and infall, the 
gas replenishment time is of order 3. 10 7 yr. The inferred duty cy- 
cle for AGN activity and nuclear star formation is a few percent. 
For simplicity, we assume the cloud to be spherical, although the 
physical mechanisms at work are expected to similarly apply for 
a wide range of geometries. 

2.2. The role of AGN in cloud disruption 

Let us consider momentum flux balance for a cloud of 
mass M c i 0Ul j and internal velocity dispersion o~. d{M c [ ou d(f)jdt = 
M w v w or Mtfoud = fEfgLEtdynico-y 1 , where a geometrical factor 
f g ~ 1 allows for the jet inefficiency in driving a quasispherical 
bow shock, Le is the Eddington luminosity, /e is the Eddington 
ratio (Lagn/Le) and v w ~ (0.1 - 0.3)c is the initial wind ve- 
locity for a wind of mass outflow rate M w . During the active 
phase of the AGN, we take /e = 0.1/£,o.i with /e,o.\ ~ 1, and 
one can disrupt some 10 6 M Q . The bow shock radius r s is set 

by Anp g r]vl lf = M w v w , or r » 2pc [f Ej ojfgM B H,4n 5 'v^] , 

where v,„/ = 50kms _1 v,„/,5o, n = 10 5 n5crrT 3 , and Mbh = 
4xlO 6 M BH , 4 M 

This is a conservative estimate, since the increased area of 
the thermal blast wave means that the jet-induced ram pressure 
plus ambient medium shocked gas thermal pressure exceeds the 
jet input momentum by one or even two orders of magnitude: eg. 
f g ~ 100 at ~ 10 6 yr, depending on jet evolution dWagner et al.l 
120121) . The swept-out radius is of order a few parsecs, compa- 
rable to the molecular cloud scale. We assume the cloud feeds 
the AGN (and forms stars) over an initial dynamical time, and 
so the cloud should be disrupted by the AGN over a time-scale 
oc r 4 , of order 10 6 yr. A more detailed jet-driven bubble model is 
described below. 

2.3. Star formation triggering by AGN 

BH feeding and triggered star formation occur simultaneously. 
There is overwhelming evidence that AGN are capable of 
quenching star formation. This is certainly likely to be the case 
for the GMC within a few parsecs of the central black hole. 
However there is also expected to be star formation within the 
self-gravitating accretion disk that directly feeds the SMBH via 
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Figure 1. Schematic model of jet/blast wave - cloud interaction 
in central 10 pc. The jet inflates a cocoon which punches into 
the GMC. The cloud is subjected to the wind arising from the 
systematic expansion of the cocoon, which has the dual effect of 
compressing its outer layers, thereby triggering star formation, 
and accelerating some the newly formed stars to become hyper- 
velocity stars. 



both fragmentation and triggering. We expect the triggered nu- 
clear star formation rate to be regulated by the enhanced pres- 
sure and to be propor tional to the square root of the pressure 
dSilk & Normanl2009l) . 

For a SFE of 0.02 per dynamical time in a 1O 6 M GMC, one 
forms 2.1O 4 M of stars per dynamical time, 10 5 yr, or 20 times 
the nuclear star formation rate (SFR). The mass fraction of this 
GMC that forms stars over a cloud lifetime is of order ~ 1%. 
The high star formation rates observed near AGN motivate us 
to assert that the star mass fraction formed is elevated, to say ~ 
10%. The enhanced SFR is justified if the pressure of the central 
cloud is elevated by a factor ~ 100 compared to nearby GMCs, 
as is the case if the cloud is self-gravitating and of size 10 pc and 
mean density ~ 10 4 ~ 5 cm~ 3 . 

A top-heavy IMF is motivated observati onally for the star- 
forming disks at the GC (iBarfko et alJl2010h . If the cloud con- 
tinues to form stars for ~ 10 cloud free-fall times, it forms some 
20% of its mass in stars over its lifetime, or 2.10 5 M Q in stars. 
There are of order 10 4 OB stars for a top-heavy IMF. 



3. The model 

The si mulations of the inte raction of a jet with a gas-rich disk 
galaxy iGaibler et al.l d201 ll) show that radio lobes develop and 
drive a quasispherical blast wave, as long as the jet length is 
smaller than the vertical scale height of the disk. The jet then 
escapes vertically and the sideways expansion of the blast wave 
stalls. HVS form in the shell at the sites of pre-existing over- 
densities, as long as the shell expands rapidly. After breakout, 
more stars form with smaller velocities. 

This picture may also apply to the disky cloud around our 
galaxy's SMBH. In this case, the jet should have a total power L 
of a fraction of the Eddington luminosity of Le = 5.10 44 ergs/s 
(for mass 4.1O 6 M ). The jets will c ollimate and develo p lobes 
from around the inner scale L\ (e.g. iRrause et al. given 
by Li ~ 0.01 pc (L/L E ) 1/2 (po/10 5 mp)- 1/2 (vj/c)- 3/2 , where the 
central cloud density is po, and the jet velocity is Vj. 

Within a molecular cloud centred on the SMBH, a jet should 
develop lobes and bow shocks just as in t he more familiar ex- 
tragalactic radio sources, and the results o f IGaibler etal]d20 111) 



may be scaled down to the GC, due to the high density of molec- 
ular clouds. A sketch of our model is provided in FigureQ] 

The ejection efficiency is the central idea of our model. With 
2% star formation efficiency (per dynamical time), and a cloud 
life-time of 10 dynamical times, regulated by AGN disruption, 
one could form a total of 10 4 stars each of mass ~ 2OM (assum- 
ing a top-heavy IMF). For a more conventional IMF, one would 
form far fewer OB stars, perhaps only ~ 100, requiring a corre- 
spondingly higher ejection efficiency. We infer that ~ 1% of OB 
stars need to be ejected in the event. In this case, the total stellar 
mass ejected is ~ 2OOOM and the associated gas mass ejected 
is 4xlO 4 M , or a few percent of the GMC. 

These numbers seem to work for the Milky Way Galaxy. 
Let us assume that in one of these events, 100 stars of aver- 
age mass 1OM are ejected at 300 km/s. Their total momen- 
tum flux is 6xl0 43 g cm/s. Over 10 6 yr, the momentum flux 
ejected is 2xl0 30 dynes. The SMBH at the GC has an Eddington 
luminosity of 5xl0 44 ergs/s (for mass 4xlO 6 M ) and we have 
assumed that it radiates at 10% of the Eddington rate. If the out- 
flow were momentum-driven, the ratio of AGN radiative pres- 
sure LEdd(4-7rr 2 cy l to mechanical wind pressure Mv(4nr 2 y l is 
LEdd(Mvc)~ l = r/c/v ~ 10, where v is the wind velocity and 
rj ~ 0.1 is the radiative efficiency of accretion. The associ- 
ated momentum flux (for optical depth unity, f=l) is fLsdd/c = 
1.7 x 10 34 /dynes (actually observatio ns require a factor / ~ 10 
to account for the M BH - cr relation dSilk & Nusseill20ldl) '). In 
this case, the ejection efficiency is 1.2 x 10~ 4 //. This seems rea- 
sonable even for fe ~ 0.1. 

The momentum transfer is determined by the ram pressure 
of the jet-driven expanding blast wave. The increased surface 
area of the expanding bow shock results in a sub stantial boost 
in momentum-driving dKrause and Gaiblerll2010l) . more quanti- 
tatively demonstrated in lWagner et al.l d2012l) . Additional ways 
of enhancin g the momen t um transfer include an energy-drive n 
blast wave dkind d2003l) : iFaucher-Giguere & Ouataertl d2012h) 
and a non-isofhermal dark halo dMcOuillin & McLaughlin! 
120121) . Similar enhancements are requi red to expla i n the ob- 
served outflows, as recently advocated in lSturm et al.l (1201 ll) . 

We estimate the ejection velocity of a gas clump by assum- 
ing that a jet with power 10 44 L44ergs/s drives a blast wave into 
a cloud with initial density 10 5 «5 cirT 3 and radius 5Rs pc (we 
ignore density stratification here). We emphasize that it is not 
actually the jet itself that accelerates clumps and induces star 
formation, but rather the associated cocoon/shock wave. Note 
also that the tidal radius is r oughly r t ~ r c (M c /Msj/) (e.g. 
dMurrav-Clav & Loebl (1201 lh l. If the tidal radius is larger then 
the extent of the cloud, star formation is not affected by tidal 
forces, and self-gravitating cores are further stabilized. 

The bubble radius evolves as R b = (5L/(47rp )) 1/5 f 3/5 , 
and the cloud is eroded over a time R 5 J 3 ((47r/5)p/L) 1 ^ ~ 



5xl0 J R 



3r,5/3, 



. 5 > 5 /L 4 4) 1/3 yr. 
star formation in the whole cloud abruptly terminates, the dif- 
fuse gas of the cloud is compressed into a dense shell that 
cools on a very short timescale t c « lQQT^ 2 n 5 l years, where 
10 7 T7K is the initial temperature of the shocked cloud gas. 
Any initial density inhomogeneitie s will then clum p due to thin 
shell and gravitational instabilities dVishniaclll983l) . At the edge 
of the cloud, the clumps should acquire a velocity of v& = 

1 /3 —1/3 —2/3 

600L 4 4 « 5 R 5 km/s. The momentum of the shell when it 
reaches the edge of the cloud would be 10 47 g cm/s, sufficient to 
explain the HVS observations, even if only a fraction of the gas 
forms stars. The time scale difference between acceleration and 



shorter than the free fall time. While 



3 



Joseph Silk et al,; Jet interactions with a giant molecular cloud in the Galactic centre and ejection of hypervelocity stars 



fragmentation is important: the cloud will in fact be disrupted, 
but stars are formed in the process . 

These clumps should cruise at basically constant speed until 
they form stars, because due to clumping and density stratifi- 
cation the clumps are much denser than their surroundings. If 
we assume a ten to hundredfold increase of the density in the 
shell, compared to the average cloud density, we infer a free fall 
time of several 10 4 up to 10 5 yrs. Assuming an exterior density 
of 10 cm 4 , as observed in the hot gas today (Ba ganoff et alj 
120031) . the expanding shell would have to accelerate to several 
1000 km/s. This triggers Rayleigh-Taylor instabilities, so that the 
clumps disconnect from the shell, which reforms at higher speed 
and continues to sweep away the interstellar medium. Even at 
the present epoch, the observed jet might lead to some HVS ejec- 
tions, as direct jet-cloud interaction s should lead to molecul ar 
cloud acceleration up to « 270 km/s (lYusef-Zadeh et alj]2012l) . 

HVS formation might also manifest an anisotropic geome- 
try. There is a suggesti on that the hyp ervelocity stars are found 
in two thin disk planes dLu et al.ll2010l) . A disk-like geometry is 
indeed sug gested by the nume rical simulations of triggered star 
formati on dGaibler et alj|201 Skewed jets a re commonly ob- 
served dKinnev et al.ll2000t lLagos et all 1201 ll) . and it is hardly 
plausible that successive mis aligned feeding disks wo uld prefer 
the same plane of symmetry (lLodato & Pringlell2006l) . 

Another manifestation of a recent AGN triggering episode 
may be the recently discovered Fermi gamma ray bubbles. The 
Fermi bubbles require E bubble = 10 55 erg or more in energy in- 
put into the relativistic plasma. The associated momentum is 
E bubble I c = 3xl0 44 g cm/s. Over 10 6 yr this amounts to a rate 
of 10 31 dynes. It is also roughly equal to the momentum flux in 
the ejected stars for the usual SFE of 10% of that in the total gas 
plus stars ejected). This seems reasonable if most of the energy 
injected is thermal. In fact the required efficiency for driving the 
star-forming clum ps is even smaller as the AGN wind model for 
the Fermi bubbles (IZubovas & Navakshinl2012l) requires ~ 10 57 
ergs injection in a few 10 5 yr and a jet power of ~ 10 44 ergs/s. 

4. Cosmological implications 

One consequence of jet-induced clump acceleration leading to 
young star ejection is that jet-induced cloud/star motions could 
retain some orbital memory of a central as opposed to a disk in- 
jection mode. Another global consequence is that nuclear ejec- 
tion would have happened more often in the past for our MWG, 
and much more vigorously in the past for galaxies with more 
massive BH than that of our MWG. 

Some 100 1OM stars ejected per 10 7 yrs imply 10 5 star ejec- 
tions over the age of the MWG. Let us allow a factor of ~ 10 
for past enhanced activity, that is some 10 6 stars or 1O 7 M are 
ejected from the MWG. These stars would become SNII in the 
intergalactic medium (IGM). The inferred mean metallicity gen- 
erated is 10~ 4 solar (since a L* galaxy of 10 n M o accounts for 
solar yields). This provides the observed IGM abundance floor. 
Moreover the IGM abundance floor (which is indeed about this 
value from Ly alpha forest data) should be ff-enhanced. The 
stars would have travelled a distance (300- 3000) km/s x 6 Myrs 
or 1.8-18 kpc. Over 3xl0 7 yrs, this amounts to 9-90 kpc. This 
results in IGM enrichment on halo and on group scales. This 
is about what might be expected, for an enrichment of the sur- 
rounding MWG IGM at z=0, and will provide ubiquitous en- 
richment for the IGM at z ~ 2 when AGN activity peaks and 
distances are smaller. 

Of course, the number of events depends on the gas refu- 
elling rate. Since the gas accretion and SMBH fuelling rates 



were higher in the early MWG, one ejects even more stars. A 
rough estimate comes from scaling the SFR, which is fed by 
cold gas accretion and also scales with the BH ac cretion rate 
by a ratio of about 1000 (see iMullanev etalJ (l2012h "). The early 
SFR is enhanced by about a factor of 10 in order to account for 
disk chemical evolution, eg. ICescutti et all d2009t) . This gives a 
factor of approximately 10 boost in OB star ejection relative to 
a constant rate, with corresponding implications for the time- 
dependence of the chemical evolution of the IGM. The fact that 
the triggering mode is more efficient at high z suggests that any 
surviving high velocity stars from early episodes should be bi- 
ased to lower metallicity. 

A time delay of 50-100 Myr has been measure d for several 
hyper velocity stars apparently ejected from the GC dBrown et al.l 
2012). These represent only a small fraction of the hypervelocity 
stars. Our model predicts that only ~ 1% of the newly formed 
jet-induced stars are ejected promptly. It is possible that a sig- 
nificant fraction of the newly formed massive stars may have 
massive binary companions, leading to a time delay of 4xl0 7 
yr or longer, if the binary companion exploded as a core col- 
lapse supernova or as a prompt SN Ia single degenerate core 
(M aoz. Mannucci a nd Brandt 1 120121) . 

In summary, we have argued that many of the observed 
hypervelocity stars were generated by the interactions of an 
AGN jet-driven cocoon from the central black hole with a dense 
molecular cloud. There are broader implications for nuclear star 
formation, the enrichment of the intergalactic medium and the 
possible origin of the Fermi bubbles. 
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